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ABSTRACT 

We investigate the stellar populations in the star forming ring of the luminous infrared galaxy 
NGC 7469. We use Hubble Space Telescope multi-wavelength (UV through NIR) imaging comple- 
mented with new if-band ground-based long-slit spectroscopy, and mid-IR and radio maps from the 
literature. Spectral energy distributions (SEDs) and evolutionary synthesis models have been used to 
characterize the star formation at different scales from those of individual star clusters (tens of pc) 
to that of the entire star- forming ring (kpc scale). At the smallest scales two different populations of 
massive (1 — 10 x 10 6 M Q ) clusters are identified. About 25% of the clusters are young (1 — 3Myr) 
and extincted {Ay ~ 3 mag), whereas the vast majority are of intermediate age (~ 9 to 20Myr) 
and less obscured (Ay ~ 1 mag). At larger (hundreds of pc) scale, an analysis of the integrated 
SED and spectroscopic data of the ring indicates the presence of two stellar populations. The young 
(5 — 6 Myr) and obscured stellar population accounts for the Br7 emission and most of the IR lumi- 
nosity, and about one-third of the stellar mass of the ring. The much less obscured intcrmcdiate-age 
population has properties similar to those of the majority of the (older) 1.1 /im-selected star clus- 
ters. The distribution of these two populations is clearly different and even spatially anti-correlated. 
The UV-optical-NIR continuum (including the majority of the clusters) of the ring traces mostly the 
mildly obscured intermediate-age population, while the MIR and radio peaks mark the location of the 
youngest and obscured star-forming regions. Moreover, the two brightest MIR and radio peaks are 
spatially coincident with the ends of the nuclear molecular gas bar. This study emphasizes the need 
for multi-wavelength, high-angular resolution observations to characterize the star formation in the 
dust-obscured regions commonly present in LIRGs. 

Subject headings: galaxies: individual (NGC 7469) — galaxies: clusters: individual — galaxies: star 
clusters — galaxies: starburst — galaxies: stellar content — infrared: galaxies 



1. INTRODUCTION 

Star formation (SF) in Luminous Infrared Galaxies 
(LIRGs, 10 11 < iiR[g_iooo/im] < 10 12 L Q , see Sanders & 
Mirabel 1996, and Lonsdale, Farrah & Smith 2006 for de- 
tailed reviews) is one of the most energetic phenomena in 
the Universe. In a large fraction of local LIRGs the bulk 
of their SF is generated within their innermost regions 
(a few kpc), in ring- like structures, mini-spiral arms, and 
regions of compact emission (Soifer et al. 2001; Alonso- 
Herrero et al. 2002, 2006a). These regions of SF contain 
large amounts of molecular gas (> 10 8 M Q , Scoville et 
al. 1991 Downes & Solomon 1998; Gao & Solomon 2004) 
and dusty, young (< 100 Myr), massive (10 5-7 M Q ) super 
star clusters (SSCs; see, among others, Surace, Sanders, 
& Evans 2000; Scoville et al. 2000; Alonso-Herrero et al. 
2000, 2002; Wilson et al. 2006). The intense SF is often 
related to past interactions suffered by the host galaxy 
with a close companion and/or to the presence of a bar 
(see Knapen 2004 for a review). Both mechanisms are 
believed to be efficient in driving large amounts of gas 
towards the nuclear regions. 

NGC 7469 is a SABa galaxy with a Seyfert 1 nucleus, 
located at a distance of ~ 65 Mpc (cz ~ 4917 km s -1 ; 
Heckman et al. 1986; H = 75 km s" 1 ). A companion 
galaxy, IC 5283, located at ~ 22 kpc (Burbidge, Bur- 
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bidge & Prendergast 1963; Dahari 1985), is believed to 
have interacted with NGC 7469 more than 1.5 x 10 8 yr 
ago (Genzel et al. 1995, hereafter G95). NGC 7469 
harbors a circumnuclcar star-forming ring with an ap- 
proximate diameter of 5" (~ 1.6 kpc) which accounts 
for approximately two-thirds of the bolometric luminos- 
ity of the galaxy. This circumnuclcar emission was first 
detected in the radio and in the near-infrared (NIR) (Ul- 
vestad, Wilson, & Sramek 1981, and Cutri et al. 1984, 
respectively). The ring was later observed in the opti- 
cal (Wilson et al. 1986; Mauder et al. 1994; Malkan, 
Gorjian, & Tarn 1998), NIR (G95), mid-infrared (MIR, 
Miles et al. 1994; Soifer et al. 2003; Gorjian et al. 2004), 
far- infrared (FIR, Papadopoulos et al. 2000), and again 
at radio wavelengths (Meixner et al. 1990; Wilson et al. 
1991). NGC 7469 has a large-scale (several kpc) stellar 
bar detected in the NIR (Knapen, Shlosman & Peletier 
2000) that gives the galaxy its 'SAB' classification, but 
also harbors a nuclear gas bar (with a size similar to the 
diameter of the ring) that has been reported by various 
authors (Laine et al. 2002; Davies, Tacconi, & Genzel 
2004). A bright radio supernova (SN2000ft) has been 
discovered in the circumnuclear ring (Colina et al. 2001; 
Alberdi et al. 2006), and recently its optical counterpart 
has been identified (Colina et al. 2007). 

G95 presented a combination of NIR broad-band 
speckle and Fabry-Perot imaging, as well as integral-field 
spectroscopy of the nuclear regions of NGC 7469 at sub- 
arcsecond resolution. They were able to separate the 
emission from the unresolved Seyfert nucleus from that 
of the surrounding starburst ring. By comparing multi- 
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wavelength diagnostics for the ring as a whole with their 
starburst model, they concluded that SF in the ring has 
either: (a) been progressing at a constant rate for the 
past several 10 7 yr, but with few high-mass stars; or (b) 
been decaying exponentially since the onset of a burst 
15 Myr ago, with a fairly 'normal' initial mass function 
(IMF). However, they were unable to assess the possible 
variation of the SF history along the ring. 

In this paper we present Hubble Space Telescope (HST) 
UV through NIR imaging, as well as low resolution (R ~ 
600) A"-band spectroscopy at four different locations 
along the circumnuclear star- forming ring of NGC 7469. 
Our goal is to study in detail the SF properties of in- 
dividual star clusters (spatial scales of tens of parsecs) 
detected in the ring (see, e.g., Scoville et al. 2000), of 
ring sections (spatial scales of a few hundred parsecs), as 
well as of the ring as a whole. This work is organized 
as follows: in § [2] we present imaging and spectroscopic 
data and their reduction procedures and analysis; in § [3] 
we describe the fitting method employed to obtain the 
parameters of the stellar populations; in § |||we present 
the results for the smallest scales (individual clusters) 
whereas results for the ring sections and the whole ring 
are given in § [5] finally, in § [5] we discuss the results and 
in § [7] briefly summarize the main conclusions. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. HST Archival Imaging 

We have retrieved HST archival images of the cen- 
tral region of NGC 7469 taken through eight filters us- 
ing three different instruments: WFPC2, NICMOS, and 
ACS. The images cover a broad wavelength range going 
from the UV (0.218 fxm) to the NIR (2.22 fim). Detailed 
information about the HST observations is given in Ta- 
ble Q] The UV and optical images were calibrated em- 
ploying the On The Fly Recalibration (OTFR) system. 
We reduced and calibrated the NICMOS images as de- 
scribed by Alonso-Herrero et al. (2000). We subtracted 
theoretical PSFs generated with the TinyTim software 4 
(Krist et al. 1998) for each filter in an attempt to min- 
imize the Seyfert nucleus contribution when analyzing 
the properties of the ring of SF. The HST images of the 
central ~ 1.6 kpc of NGC 7469 are presented in Fig. |TJ 

2.1.1. Star Cluster selection 

Star clusters were selected at 1.1 ^m (NICMOS F110W 
filter) because of two main reasons: (a) the NIR spec- 
tral region is less affected by dust extinction than UV 
and optical wavelengths: A uv (0.22/™) ~ 3 x Ay ~ 
20 x A 1AfJ , m and A B{0 ^ im) ~ 2 x A v ~ 13 x 
Ax.i Mm (Calzetti 1997, CAL hereafter); (b) the F110W 
filter has the best spatial resolution among the NICMOS 
filters (see Table HJ. 

The positions of the star clusters in all the HST im- 
ages were determined using a procedure similar to the 
growing region method employed by the HiiPhot soft- 
ware (Thilker, Braun, & Walterbos 2000). This search 
method is also similar to that used by Sextractor 
(Bertin & Arnouts 1996). A 2a threshold above a back- 
ground selected in four regions in the outskirts of the 
ring was imposed for cluster identification. However, 

4 http: / /www. stsci . edu/sof tware/tinytim 



this background was in some cases not representative 
of the strong, complex and variable "background" emis- 
sion surrounding the clusters (that is, underlying unre- 
solved emission within the ring acting as a pedestal back- 
ground) . The search method had to be stopped manually 
when noise peaks began to be identified as positive de- 
tections in the images, specially in those with the lowest 
S/N ratio (F547M and F814W). A total of 30 clusters in 
the ring of SF was selected at 1.1 /xm. Their locations are 
given in Tabled] (see also Fig. |2J). Out of the 30 1.1 /xm 
star clusters within the ring selected by our method, only 
five (CI, C3, C5, C6, and C7) were identified by Scoville 
et al. (2000) using the same NIR dataset. 

In order to construct spectral energy distributions 
(SEDs) of the 1.1 /xm-selected star clusters, we cross- 
correlated their positions with those of clusters detected 
at other wavelengths. Clusters located at distances < 
1.5 x FWHMpnow (< 0'.'175 ~ 55 pc) were identified as 
counterparts of the 1.1 /xm clusters. Because the F330W 
and F606W images have both good S /N and higher spa- 
tial resolution than the F110W image, a few 1.1 /xm clus- 
ters were found to have multiple cross-identifications at 
these wavelengths. Table [3] provides information about 
cluster identification and cross-correlation. For all the 
images except those with the highest spatial resolutions 
(F330W and F606W), most of the clusters identified at a 
given wavelength are indeed those also selected at 1.1 /xm, 
even in the UV. The fraction of the F330W and F606W 
clusters cross-identified at 1.1 /im is smaller than in the 
other filters, but consistent with the larger number of 
clusters detected in these two filters. 

2.1.2. Cluster photometry 

The cluster intensity was modelled with a 2D Gaussian 
function whose width was fixed to that measured from 
theoretical PSFs generated with TinyTim (Krist et al. 
1998) for each filter. Additionally, our method fits the 
emission surrounding the source with a 2D plane, which 
is then interpolated to the region where the cluster is lo- 
cated. The region used to fit the background was defined 
as the area enclosed by a two pixel wide ring with its in- 
ner radius containing the 90% of the TinyTim PSF flux. 
This is done after rejecting those pixels that exceed the 
fitted background by 2.75<r (i.e., pixels in the background 
with possible contamination from other clusters). A \ 2 ~ 
minimizing Marquardt method (Bevington & Robinson 
2003) implemented in IDL was used to derive accurate 
positions and fluxes of the point sources. However, since 
the -ffSTPSFs are not Gaussians, an aperture correction 
(filter-dependent) was needed. These corrections were 
obtained by modelling the individual theoretical PSFs 
generated with TinyTim in the same way as the star 
clusters. 

In order to estimate the photometric uncertainty of this 
method, we generated a grid of test images where Tiny- 
Tim PSFs with random positions and variable intensi- 
ties were placed. The number density of point sources 
was approximately 0.01 pixel -2 . PSFs separated by 
< 2 — 3 pixels were common in the test images, sim- 
ilar to the the cluster density observed in the ring of 
NGC 7469. We ran simulations varying the background 
intensity and standard deviation, as well as the PSF in- 
tensity with respect to that of the background. We ap- 
plied our fitting method to measure the PSF fluxes and 
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TABLE 1 

HST Imaging data 



Filter 


Instrument 


Detector 


Plate scale 


Date 


t (a) 


Orientation 


TT FWHM< b ) 


TT FWHM< b ) 


ID 


Pl(c) 








(arcsec pixel -1 ) 


(yy/mm/dd) 


(s) 


(degrees) 


(arcsec) 


(PC) 






F218W 


WFPC2 


WF3 


0.1 


99/06/28 


1200 


107.54 


0.176 


55 


6358 


Colina 


F330W 


ACS 


HRC 


0.026 


02/11/20 


1140 


-107.75 


0.051 


16 


9379 : 


Schmitt 


F547M 


WFPC2 


PC 


0.0455 


00/05/13 


12 


-61.36 


0.076 


24 


8240 


Wilson 


F606W 


WFPC2 


PC 


0.0455 


94/06/10 


500 


-68.35 


0.073 


23 


5479 ■ 


Malkan 


F814W 


WFPC2 


PC 


0.0455 


00/05/13 


7 


-61.36 


0.085 


27 


8240 


Wilson 


F110W 


NICMOS 


NIC2 


0.0755 


97/11/10 


352 


-150.15 


0.117 


37 


7219 : 


Scovillc 


F160W 


NICMOS 


NIC2 


0.0755 


97/11/10 


352 


-150.15 


0.154 


48 


7219 : 


Scovillc 


F222M 


NICMOS 


NIC2 


0.0755 


97/11/10 


480 


-150.15 


0.202 


64 


7219 : 


Scovillc 



Note. — (a) Total exposure time of the final image; (b) FWHM of TinyTim PSFs; (c) Proposal ID and PI. 



TABLE 2 

Positions and flux densities of the 1.1 /xm-selected star clusters in the ring of SF 



# (a) 




y(b) 


F218W( C ) 


F330W( C > 




F547M< C ) 


F606W< C ) 




F814W( C ) 


FllOWW 


F160W< C ) 


F222M( C ) 


CI 


-1.54 


-0.67 


9.88x10" 


16 


1# 3.44xl0" 


16 


1.03x10" 


16 


7.73x10" 


17 


5.90x10" 


17 


3.58x10" 


17 


2.79x10" 


17 


1.93x10" 


17 


C2 


-0.96 


0.99 


2.59x10" 


16 


1.69x10" 


Ifi 


1.05x10" 


16 


8.85x10" 


17 


4.95x10" 


17 


2.63x10" 


17 


1.48x10" 


17 


9.85x10" 


18 


C3 


0.93 


1.34 


9.41x10" 


17 


7.38x10" 


17 


6.74x10" 


17 


4.94x10" 


17 


2.36x10" 


17 


2.28x10" 


17 


1.79x10" 


17 


1.22x10" 


17 


C4 


-0.74 


-1.20 






6.87x10" 


18 






i # 2.96xl0" 


17 






1.97x10" 


17 


1.93x10" 


17 


1.60x10" 


17 


C5 


-1.73 


0.10 


9.23x10" 


17 


6.42x10" 


17 


4.03x10" 


17 


4.67x10" 


17 


2.27x10" 


17 


1.80x10" 


17 


1.24x10" 


17 


7.93x10" 


18 


C6 


1.19 


1.11 


4.95x10" 


17 


i # 2. 69x10" 


17 






2.03x10" 


17 


1.63x10" 


17 


1.68x10" 


17 


1.63x10" 


17 


1.43x10" 


17 


C7 


-0.46 


-1.54 


1.29x10" 


16 


i # l. 31x10" 


16 


7.03x10" 


17 


7.59x10" 


17 


3.74x10" 


17 


1.41x10" 


17 










C8 


-0.37 


1.49 






i # 8.66xl0" 


18 






1.47x10" 


17 






1.34x10" 


17 


1.02x10" 


17 


9.85x10" 


18 


C9 


-1.28 


0.59 






2# 3.23xl0" 


17 


1.92x10" 


17 


1# 2.02xl0" 


17 






1.22x10" 


17 


1.18x10" 


17 


9.39x10" 


18 


C10 


-0.58 


1.58 






4.83x10" 


18 






1.31x10" 


17 






1.20x10" 


17 


1.23x10" 


17 






Cll 


1.27 


-0.68 






2.89x10" 


18 














1.09x10" 


17 


1.06x10" 


17 


7.67x10" 


18 


C12 


1.27 


-0.90 






2.89x10" 


18 






8.71x10" 


18 






1.06x10" 


17 










C13 


-1.79 


-0.37 






1.36x10" 


17 






1.47x10" 


17 






1.02x10" 


17 






4.35x10" 


18 


C14 


-1.11 


0.35 


1.19x10" 


16 


2# 6.62xl0" 


17 


3.16x10" 


17 


i # 3. 70x10" 


17 


1.65x10" 


17 


1.00x10" 


17 


6.88x10" 


18 


6.55x10" 


18 


C15 


-0.28 


-1.28 


3.57x10" 


17 




2.14x10" 


17 


i # 3.64xl0" 


17 


1.79x10" 


17 


8.91x10" 


18 


8.28x10" 


18 


7.05x10" 


18 


C16 


1.72 


-0.02 


3.67x10" 


17 


i # 3. 11x10" 


17 






1.75x10" 


17 






8.89x10" 


18 


:6.81xl0" 


18 






C17 


1.04 


0.37 






1.13x10" 


17 






1.17x10" 


17 






7.82x10" 


18 


6.99x10" 


18 






C18 


0.17 


-0.92 


4.22x10" 


17 


2.03x10" 


17 






2.29x10" 


17 






7.73x10" 


18 


8.36x10" 


18 


6.28x10" 


18 


C19 


1.36 


-0.39 


1.14x10" 


17 


















7.45x10" 


18 


7.56x10" 


18 


5.93x10" 


18 


C20 


-1.49 


-1.05 






















7.35x10" 


18 


8.65x10" 


18 


9.91x10" 


18 


C21 


1.52 


0.05 






3.06x10" 


18 














7.20x10" 


18 


:6.81xl0" 


18 


5.31x10" 


18 


C22 


0.91 


-0.15 


4.77x10" 


17 


3.15x10" 


17 


3.42x10" 


17 


2.59x10" 


17 


1.77x10" 


17 


6.94x10" 


18 










C23 


-1.34 


-1.37 






i # l. 18x10" 


17 






7.69x10" 


18 






6.85x10" 


18 


5.24x10" 


18 






C24 


-1.55 


0.31 


4.15x10" 


17 


l # :3.54xl0" 


17 






1.60x10" 


17 






6.18x10" 


18 










C25 


-1.34 


0.13 






i # :l. 37x10" 


17 






7.78x10" 


18 






5.97x10" 


18 


4.34x10" 


18 






C26 


0.90 


0.24 




















4.43x10" 


18 










C27 


1.29 


0.51 






2.31x10" 


18 














4.26x10" 


18 










C28 


0.19 


1.23 






4.38x10" 


18 






8.31x10" 


18 






4.09x10" 


18 


2.87x10" 


18 


3.41x10" 


18 


C29 


0.50 


0.98 






1.48x10" 


17 






1.39x10" 


17 






2.36x10" 


18 










C30 


-0.98 


-0.28 






















2.18x10" 


18 











NOTE. — (a) Cluster label. The clusters are sorted by their F110W flux; (b) Positions of the clusters relative to the nucleus (in arcsec). A positive 
sign indicates west and north directions; (c) Flux densities (in erg s _1 cm" 2 A -1 ) at each filter. The number with a hash symbol next to a flux density 
indicates the number of additional counterparts that have been detected aside from the first one. The quoted flux is the sum of all the counterparts. A 
colon next to a flux density indicates that the same cluster is associated with two different 1.1 /im star clusters. 



Resolving the Stellar Populations in NGC 7469 



5 




-2-10 1 2 

arcsec 

FlG. 2. — The contours, shown in a logarithmic scale, represent 
the HST/TSSICMOS F110W emission. The 1.1 /im-selected star 
clusters are marked with small open circles with sizes equivalent 
to the FWHM of the F110W PSF (~ 37 pc, see Table Q). The 
location and P. A. of the UKIRT/CGS4 slits are also marked, wit 
the 0'.'6 slit width and extraction lengths (see Table |4j shown to 
scale. The ring area defined by G95 (enclosed by two open circles 
of radii Of.' 7 and 2"5, see text) is also shown. North is up, and east 
to the left. [See the electronic edition of the Journal for a color 
version of this figure.] 



TABLE 3 

Photometric uncertainties, number of clusters and 
cross-identifications 



Filter 


Phot 


No. 


No. 


Fraction of 




uncert/ a ) 


clusters^ 3 ' 


cross-ID^ ' 


1.1 (im selected'" 1 ' 


F218W 


10% 


16 


13 


81% 


F330W 


5% 


56 


25+11 


45+20% 


F547M 


15% 


10 


9 


90% 


F606W 


5% 


12 


23+4 


55+10% 


F814W 


12% 


11 


10 


91% 


F110W 


8% 


30 


30 


100% 


F160W 


8% 


20 


20 


100% 


F222M 


9% 


18 


17 


94% 



NOTE. — (a) Average estimated uncertainty for the cluster photometry 
(see text for details); (b) Number of identified clusters; (c) Number of 
cross-identified 1.1 /im clusters in each image plus the number of clusters 
identified as double or triple counterparts for the F330W and F606W 
filters; (d) Fraction of cross-identified 1.1 /mi clusters with respect to the 
total number of clusters detected in each image. 



compare them with the input values. Our point source 
modelling method was compared with two other classical 
photometric tools: aperture photometry, and the tasks 
of daophot (Stetson 1987) in iraf 5 . We found that our 
method is capable of recovering fluxes with an accuracy 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



TABLE 4 
UKIRT/CGS4 slit parameters 



Slit 


PA. 


f (a) 
t-exp 


Aperture( b > 


Aperture( b ) 




(degree) 


00 


(arcsec) 


(PC) 


Al 


123.9 


5280 


0'.'61 X 2'.'40 


192 x 756 


A2 


107.7 


5520 


0'.'61 x 2^93 


192 x 923 


A3 


42.1 


6240 


0'.'61 x 3'.'00 


192 x 945 


A4 


174.7 


5280 


0'.'61 x 2^69 


192 x 848 


Nucleus 




960 


0'.'61 x l'.'OO 


192 x 315 


Note. — 


(a) Exposure time; 


(b) Dimensions 


of extraction 



apertures. 



typically as high as that of daophot and in some cases 
twice as high (depending on the filter and the parame- 
ters of the simulation), and ~ 2 — 4 times higher than the 
classical aperture photometry method. This is probably 
due to the detailed background modelling performed by 
our method. 

These simulations were also used to estimate the un- 
certainties (see TableEJ of the cluster photometry in each 
image by selecting the simulation that best matched the 
background parameters (standard deviation and relative 
intensity of the clusters with respect to the background) 
of the ring. The photometry for the 1.1 /im-selected star 
clusters is given in Table [2l In the case of multiple cross- 
identifications, the sum of the flux densities of all cross- 
identified clusters at a given wavelength was used as the 
representative of the cluster selected at 1.1 /im. 

2.2. Long-slit NIR spectroscopy 

Long-slit spectroscopy of the central region of 
NGC 7469 was obtained with the common-user NIR 
spectrograph CGS4 (Mountain et al. 1990) on UKIRT 
during 1999 October 5 and 6. We employed the 
0'.'61 pixel^ 1 plate scale, a 1 pixel wide slit, and the 
40 line mm -1 grating in first order, which delivers com- 
plete spectral coverage from 1.85 — 2.45 /im at a resolving 
power of 600. Conditions were photometric with 0'.'6 see- 
ing throughout the first night, but increasing amounts of 
cirrus curtailed observations on the second night. A total 
of 4 slit position angles and nuclear offsets were selected 
(see Fig. [2]). Details of the observing run are given in 
Table 2J After setting the slit to the desired position an- 
gle, a manual search for the if-band nuclear peak was 
conducted, following which accurate offsets were applied 
by use of the UKIRT crosshead to place the slit right at 
the desired locations. In addition, a single slit position 
passing right through the nucleus was observed. 

The data reduction was accomplished using a combi- 
nation of tasks within the Starlink CGS4DR and NO AO 
iraf packages, as described in Ryder et al. (2001). The 
relatively compact nature of the infrared-emitting central 
region of NGC 7469 allowed sky subtraction to be per- 
formed by sliding the region of interest a total of 38'.'5 arc- 
sec along the slit in between 'object' and 'sky' exposures. 
The difference of each object-sky pair on a given slit posi- 
tion was co- added into groups. The extraction apertures 
were defined independently for each grouped image, in 
consultation with Fig. [2l Owing to the still somewhat 
coarse spatial resolution of CGS4 when compared with 
the typical size and separation of the clusters in the star- 
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forming ring, it proved impossible to reliably extract the 
emission from individual clusters. Instead, the emission 
from groups of clusters has been summed, representing 
the (luminosity- weighted) mean emission from an ensem- 
ble of clusters. This nevertheless allows us to compare 
the emission, and thus the SF history, from four distinct 
regions along the circumnuclear ring. 

2.2.1. Photometric calibration 

Since the weather conditions were not photometric 
during the second night, we made use of the available 
HST/NICMOS F222M (~ K) image to perform the pho- 
tometric calibration of the UKIRT spectra. A TinyTim 
PSF scaled to best match the bright nuclear emission was 
subtracted in order to minimize its contamination to the 
circumnuclear region. We simulated the sizes and orien- 
tations of the UKIRT extraction apertures (Table [4]) on 
the F222M image (smoothed to the spatial resolution of 
the spectra) and measured the flux densities within the 
four slits. In a similar fashion we measured the fluxes of 
the spectroscopic apertures on the rest of the HST im- 
ages in order to obtain the SED of each ring section. The 
uncertainties of flux densities were estimated by shifting 
the positions of the simulated apertures ± 1 pixel in the 
x and y directions over the HST images, and also by 
taking into account flux variations within the simulated 
apertures when we manually chose various intensities to 
scale the TinyTim PSF with the nuclear peak. Flux den- 
sities of the spectroscopic apertures, together with their 
uncertainties, are given in Table [5] We also obtained 
integrated photometry through all the HST filters (see 
Table [5]) of the circumnuclear ring region of NGC 7469 
for the area enclosed by two circular apertures of radii 
ft' 7 and 2'.'5, following the analysis of the ring done by 
G95. 

2.2.2. Spectroscopic analysis 

Fig. [3] shows the spectra of the combined data for each 
of the four ring sections (apertures), together with the 
nuclear spectrum. The main spectral features present 
in all spectra are marked at the top. Line fluxes at 
A < 2 /im have larger intrinsic uncertainties due to at- 
tempts to correct for atmospheric transmission, which 
decreases about 20-40% from 1.8 — 1.96 /xm (these lines 
are marked with an asterisk in Table[6]). Although no sys- 
tematic under/over-correction of the continuum seems to 
have been introduced, these values should be taken with 
caution. 

In order to define a consistent continuum right out into 
the region of extensive CO absorption beyond 2.3 /im, we 
fitted a power-law of the form F\ cx X 13 to featureless sec- 
tions of the spectrum near 2.09 and 2.25 /im (rest wave- 
length). After normalizing each spectrum by this fit, the 
equivalent widths (EWs) and fluxes of the most promi- 
nent emission features were measured, together with the 
spectroscopic CO sp index, as defined by Doyon, Joseph 
& Wright (1994), and the photometric CO index of Fro- 

gel et al. (1978): CO p h = — "[as • Line nuxes an d 
EWs, as well as the CO indices are presented in Table [51 
The four spectra (see Fig. [3]) taken at different posi- 
tions along the circumnuclear ring are all rather similar 
both in terms of line fluxes and their EWs. The nuclear 
spectrum differs from that of the circumnuclear ring pri- 




2 2.2 2.4 



Wavelength (yttm) 

Fig. 3. — Spectra of the four UKIRT slit positions along the 
circumnuclear ring of NGC 7469, as well as through the Seyfert 
nucleus. The most prominent emission features are labelled, as 
well as the CO band. The peak at 1.945 fim in each spectrum 
(marked with an asterisk) is spurious, and arises from Br<5 absorp- 
tion intrinsic to the standard star BS 8738. 



marily in two respects: first, the extremely red contin- 
uum slope; and second, the broad hydrogen recombina- 
tion lines. We can use these distinguishing characteris- 
tics to place an upper limit on the extent of contamina- 
tion from the bright nuclear source on the circumnuclear 
spectra. A contribution of 2% or more to the integrated 
.fT-band flux from the nuclear light would be easily no- 
ticeable and would have significantly steepened the flat 
continua of the ring spectra. The absence of even such a 
small contribution means that the 0'.'6 slit and UKIRT's 
active optics facility have done an excellent job at spa- 
tially resolving the circumnuclear starburst ring from the 
Seyfert nucleus, despite the latter being only 1'.'5 from, 
but ~ 3 mag brighter than, the region sampled by each 
slit position. 

3. EVOLUTIONARY SYNTHESIS MODELS 

3.1. Modelling parameters 

We used Starburst99 (SB99 hereafter; Leitherer et al. 
1999; Vazquez & Leitherer 2005; with nebular compo- 
nent included) to model the SF properties of the 1.1 ^tm- 
selected star clusters, the ring sections, and the ring as 
a whole. The main input parameters are: type of SF, 
initial mass function (IMF), stellar evolutionary tracks, 
and metallicity, which was set to solar (a reasonable as- 
sumption for nuclear rings, Mazzuca et al. 2006; Allard 
et al. 2007). As the SF processes taking place in the 
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TABLE 5 

Observed flux densities of the spectroscopic apertures and the whole ring 



Region 


F218WM 


F330W< a ' 


F547M( a ) 


F606W< a ' 


F814W< a ) 


F110W< a ) 


F160W< a ) 


F222M( a ) 


Al 


16.5±1.0 


13.7±0.5 


10.5±0.2 


11.7±0.2 


8.6±0.3 


6.9±0.3 


5.5±0.2 


3.3±0.1 


A2 


3.4±0.3 


4.9±0.3 


6.2±0.2 


7.5±0.2 


6.4±0.2 


6.2±0.2 


5.4±0.2 


3.3±0.1 


A3 


11.6±0.8 


11.8±0.5 


10.8±0.3 


11.6±0.3 


9.61±0.2 


7.1±0.4 


5.8±0.2 


3.4±0.1 


A4 


4.5±0.4 


5.8±0.5 


6.8±0.3 


8.9±0.4 


7.4±0.2 


7.1±0.4 


6.1±0.2 


3.6±0.1 


Ring 


60.6±2.4 


62.2±1.2 


67.9±2.0 


73.6±2.2 


63. Oil. 9 


54.3±1.6 


46.5±1.9 


27.9±0.8 



Note. — (a) Flux densities and uncertainties (10 16 erg s 1 cm 2 A 1 ) for each HST filter. 



TABLE 6 

Observed line fluxes, EWs and CO indices of UKIRT spectra. 



LineM A ( b ) Al Al A2 A2 A3 A3 A4 A4 

FluxM EW< d ) Flux( c ) EW< d ) Flux< c ) EW< d ) Flux( c ) EW( d > 



Paa* 1.8751 7.97±0.41 23.7±1.6 9.57±0.35 28.6±1.4 12.6±1.0 36.6±3.0 9.94±0.90 25.7±2.8 

Br<5* 1.9451 1.34±0.20 4.10±0.70 1.29±0.30 4.00±1.00 1.73±0.28 4.90±0.90 1.56±0.24 4.30±0.70 

H 2 1-0 S(3)* 1.9576 1.84±0.16 5.70±0.50 1.91±0.24 6.00±0.80 2.25±0.23 6.70±0.70 2.33±0.26 6.50±0.70 

[SiVI]* 1.9645 0.35±0.05 1.10±0.15 0.90±0.13 2.90±0.55 0.60±0.08 1.80±0.25 0.67±0.14 1.90±0.40 

H 2 1-0 S(2) 2.0338 0.66±0.07 2.00±0.20 0.73±0.14 2.30±0.50 0.56±0.07 1.70±0.20 0.95±0.15 2.70±0.40 

Hel 2.0581 0.95±0.10 2.90±0.30 0.98±0.08 3.10±0.30 1.04±0.06 3.10±0.20 0.66±0.14 1.80±0.40 

H 2 1-0 S(l) 2.1218 1.30±0.09 3.90±0.30 1.47±0.09 4.50±0.30 1.38±0.06 4.10±0.20 1.44±0.13 4.10±0.20 

Br 7 2.1661 3.51±0.13 10.7±0.40 3.72±0.12 11.5±0.40 3.60±0.11 10.6±0.40 3.00±0.11 8.40±0.35 

H 2 1-0 S(0) 2.2233 0.15±0.05 1.45±0.15 0.47±0.12 1.45±0.40 0.39±0.07 1.16±0.20 0.36±0.06 1.02±0.20 

~CU^, OS ~ 0.21±0.05 ~ 0.23±0.05 ~ 0.23±0.05 ~ 0.20±0.05 

CO Dh 2.35 ... 0.16±0.05 ... 0.17±0.05 ... 0.17±0.05 ... 0.15±0.05 



Note. — (a) Emission lines and CO indices; (b) Rest line/band wavelength (/jm); (c) Flux (10~ 15 org s _1 cm -2 ) within each 
aperture; (d) EW in the case of emission lines (A); photometric and spectroscopic CO indices in the case of CO bands (sec text); * 
Line fluxes with large uncertainties due to poor atmospheric transmission shortward of 2/im (see text). 



nuclear regions of LIRGs and ULIRGs are believed to 
last only a few tens of Myr (Kennicutt 1998, Allard et 
al. 2006, 2007), we used instantaneous SF bursts. For 
the IMF we have used a Salpeter (1955) IMF with lower 
and upper mass limits of 0.1 M and 100 M , respec- 
tively. Also, because the star clusters located in the cir- 
cumnuclear ring of NGC 7469 are expected to be young 
(~ 1 — 50 Myr — since about 80% of the 1.1 ^im-selected 
clusters are UV-emitters — , Table [3J see also G95), we 
have used the tracks from the Geneva group because they 
better reproduce the SEDs of young, ionizing stellar pop- 
ulations (Vazquez & Leitherer 2005). 

We have constructed a grid of theoretical SEDs with 
ages ranging from 0.5 to 500 Myr. The steps of this grid 
are 0.5 Myr for the 1 — 10 Myr age interval, 1 Myr for 
10 - 20 Myr, 5 Myr for 20 - 100 Myr, and 50 Myr for 
100 — 500 Myr. The model spectra were reddened as- 
suming the extinction law of Calzetti et al. 1997 (CAL) 
in a simple dust screen configuration which is a good rep- 
resentation of the dust geometry of individual clusters. 
This dust screen model, however, may not be appropri- 
ate for extended regions where more complicated dust 
geometries are expected (see discussion by G95). 

3.2. Fitting Method 

We have compared observational and theoretical SEDs 
through a % 2 -fitting procedure (Bik et al. 2003): 

2 a \ (Jobs - /model) \ 



where N are the filters or observables in general (photo- 
metric flux densities, flux and/or EW of emission lines, 
etc.) available for each cluster, ring section, or entire 
ring, / bs and /model are the observed and model quanti- 
ties respectively, and <T bs are the weights for the fit. In 
our case, cr bs accounts for the photometric calibration 
uncertainties, as well as the uncertainties associated with 
the theoretical models. The x 2 minimization method for 
fitting SEDs produces more satisfactory results for deter- 
mining the ages, extinctions, and masses of star clusters 
detected in galaxies (see, e.g., Maoz et al. 2001; Bastian 
et al. 2005; de Grijs et al. 2005) than other widely used 
methods (e.g., color-color diagrams). 

When fitting the photometric SEDs (see next sections) 
we used flux densities instead of magnitudes to avoid 
the uncertainties associated with transforming fluxes 
based on the HST photometric system to the standard 
Johnson-Cousins (UBVIJHK) system (see discussion 
by de Grijs et al. 2005). To compare the observed flux 
densities with the outputs of the evolutionary models, 
the model spectra (which include the nebular contribu- 
tion and the main hydrogen recombination lines) were 
convolved with the eight HST filter bandpasses that take 
into account both the filter response and the camera- 
telescope system throughput. For the model uncertain- 
ties we take 5% for the optical and NIR, and 10% for the 
UV (as assumed by de Grijs et al. 2005). 

In the case of a single stellar population, there is a total 
of three parameters to be fitted: age (t), extinction (Ay), 
and stellar mass (m*). To be precise, the stellar mass is 
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TABLE 7 

Results of the x 2 -fitting of the 1.1 /iiii-selected cluster SEDs. 



Cluster 


t («) 




m Jc) x 2(d) v (e) 




(Myr) 


(mag) 


(x 10 6 M ) 



5 
5 
5 
2 
5 
4 
3 
2 
3 
1 
1 
1 
5 
4 
2 
1 
3 
1 
1 
3 
1 
1 
1 
2 



Note. — Results of the fits using SB99 models with an instantaneous 
SF burst, Geneva tracks, and CAL extinction law; (a) Age; (b) Visual 
extinction; (c) Stellar mass; (d) Best (lowest) x 2 value (see EquationfTl; 
(e) Degrees of freedom. The number with a hash symbol next to a 
cluster label indicates the number of total additional counterparts that 
have been detected among all of the images (i.e., taken into account all 
wavelengths; see Table [2] for details). 



a semi-free parameter because the fitting method sim- 
ply scales the theoretical SEDs with the mass until the 
maximum likelihood is reached. However, since the mass 
affects the solution and accuracy of the fit, it is necessary 
to include it in the equation that describes the method. 
The expected (minimum) \ 2 value of the best fit should 
be equal to the number of degrees of freedom (v = N—3), 
but in general our fits provide a higher value of x 2 (see 
Table [3 and next section). Thus, to estimate the un- 
certainty of the fitted parameters, we define v' = x min , 
where x m in ^ s the minimum value obtained for a given 
cluster. To determine the range of acceptable solutions 
we take all those solutions within Xmin ^ ^Xmim w i tn 
^Xmin = v2V. This would be equivalent to taking 
the ±lcr solutions. Although this expression is not for- 
mally correct, it provides an estimate of the acceptable 
solutions when fitting the SEDs. 

4. SMALL-SCALE STAR FORMATION: 1.1 /iM-SELECTED 
MASSIVE STAR CLUSTERS 

We have fitted the photometric SEDs of the 1.1 (im- 
selected star clusters using a single stellar population 
to derive their ages, extinctions, and masses. Since the 
number of degrees of freedom is v = N — 3, we have only 
fitted those 1.1 ^m-selected star clusters cross-identified 
in another three wavelengths or more; this resulted in a 
total of 24 fitted clusters (see Table [7]). Fig. g] shows the 
observational data and the best SED fits of the 10 most 



luminous clusters sorted by their 1.1 /jm brightness, as 
well as four fainter clusters with 3 or more degrees of free- 
dom. We also show in this figure the acceptable range 
of solutions defined as the area between the SEDs of the 
youngest and oldest valid models (within the 1 a uncer- 
tainty of the best % 2 value, as explained in § 3.2). 

The majority (75%) of clusters have SEDs best fitted 
with ages ranging between 8.5 and 20 Myr. The upper 
limits of acceptable ages for these clusters are less con- 
strained because the model SEDs in the ~ 10 — 30 Myr 
range do not vary as quickly as for younger (~ 1 — 10 Myr) 
stellar populations. We also note that clusters not de- 
tected in the UV tend to show large uncertainties in their 
fitted parameters. Nearly all of the remaining (~ 21%) 
clusters are fitted with a younger, 4=1 — 3 Myr, stel- 
lar population. Only one cluster appears to be relatively 
old (C22, 200 Myr), although with a wide range of ac- 
ceptable ages (see Table [7]). The ages of the 1.1 /zm se- 
lected clusters are similar to those of clusters detected 
in other star-forming circumnuclear rings (Maoz et al. 
2001; Alonso-Herrero et al. 2001a), although older clus- 
ters like C22 (> 100 Myr) have also been detected in 
ringed galaxies (Buta et al. 2000). 

Fig. [S] shows the age distribution of the clusters as a 
function of their position angle within the ring of SF. 
As well as in age, the two groups of star clusters also 
appear to be well differentiated in terms of their extinc- 
tion (see Table [7]). The youngest (2.5 Myr on average) 
have a mean extinction of Ay ~ 3 mag, whereas those 
of intermediate age (on average 14 Myr) tend to be less 
obscured (mean extinction of Ay ~ 1.25 mag). Some 
clusters (CI, C2, C3, C6, C8, C9, C14, C15, C18 and 
C19) present a 2.2 /xm excess (~ 20%) with respect to 
their fitted SEDs. Maoz et al. (2001) found a similar be- 
havior in clusters detected in NGC 1512 and NGC 5248, 
and interpreted it as an additional emission component 
from dust heated by the star clusters. The spatial dis- 
tribution of the 1.1 /zm-selected star clusters and their 
relation to the gas and dust distribution in the ring will 
be discussed in detail in § 6. 

A large fraction of clusters in the ring of NGC 7469 
has stellar masses between 0.5 and 10 x 10 6 M Q , al- 
though about one-third appear to be more massive (> 
10 x 10 6 Mq). These values are also in agreement with 
stellar masses inferred for SSCs found in other galaxies 
(10 4 M < M < 10 6 M for M82, de Grijs, Parmentier k, 
Lamers 2005) and especially in LIRGs (from ~ 10 4 M Q 
to few 10 7 M©; Antennae galaxy, Zhang & Fall 1999; 
NGC 3256 and others, Alonso-Herrero et al. 2000, 2001b, 
2002; Arp 299, Garcfa-Marm et al. 2006; Arp 220, 
Wilson et al. 2006). Conversely, other studies suggest 
that not all circumnuclear rings show such massive SSCs 
(Maoz et al. 1996; Maoz et al. 2001), although these 
galaxies are not classified as LIRGs. One possibility is 
that some of the 1.1 /xm-selected clusters identified as 
point sources in the ring of NGC 7469 are aggregates 
of less massive clusters, as suggested by the ACS F330W 
image (see Fig. []}. 

5. STAR FORMATION ON THE KILOPARSEC SCALE: THE 
CIRCUMNUCLEAR RING 

5.1. Stellar Populations 

The four UKIRT slits (see Fig. [5]) cover regions with 
sizes of the order of 0.2 kpc x 1.4 kpc each, increasing 



Cli# 
C2 
C3 
C4i # 
C5 
C6 1# 
C7i # 
C8 i# 
C9 3 # 
C10 
Cll 
C13 
C14 3# 
C15i # 
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Fig. 4. — Fits to the 1 . 1 /im-selected star cluster photometric SEDs. The model SEDs (SB99) use instantaneous star formation, the 
Geneva tracks, and the CAL extinction law. The ten brightest clusters are shown as well as the remaining four clusters with more than 
3 degrees of freedom (see Table are shown as (red) open circles with their 1 a uncertainties. The (blue) solid line is the best fit of the 
observational datapoints, and the open, (blue) triangles represent the photometric points resulting from convolving the best fitted SED 
with the ffi>T-filter throughputs. The (green) shaded regions represent the SED range of the youngest and oldest acceptable solutions 
(Table 0. The emission lines seen in the youngest SEDs are hydrogen recombination lines included in the modelling. As in Table [7] the 
number with a hash symbol next to a cluster label indicates the number of total additional counterparts that have been detected among 
all of the image. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 4.— 



Continued. [See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 5. — Ages of the 1.1 /xm-selected star clusters as a function 
of their P.A. within the circumnuclear ring. The position angle is 
measured so that north is degrees and east is 90 degrees. The 
shaded regions indicate the locations of the MIR/radio (see Fig. [9} 
peaks tracing the youngest stars 6 — 8 Myr) in the SF ring. [See 
the electronic edition of the Journal for a color version of this 
figure.] 



the likelihood to probe stellar populations with different 
ages. As discussed in § 4, a large fraction of the 1.1 /zm- 
selected star clusters have an average age of 14 Myr. Nev- 
ertheless, the intense B1-7 hydrogen recombination line 
emission measured from the UKIRT spectra has to be 
produced by a younger ionizing stellar population. Thus, 
to fit the photometric SEDs and spectroscopic data (Br7 
flux and its EW) of the ring sections we used a sim- 
ple combination of two stellar populations 6 . Because 
the ionizing and the MIR emitting regions are spatially 
coincident (Soifer et al. 2003 and G95, and discussion 
in § 6), the young stellar population has to be at least 
moderately obscured. Thus, we introduced an additional 
constraint for the extinction of the young stellar popu- 
lation: Ay > 3 mag, based on the typical extinctions 
over the Paa emitting regions found for local LIRGs 
(Alonso-Hcrrero et al. 2006a). From gas and stellar dy- 
namics, G95 inferred a dynamical mass for the ring of 
4.5 x 10 9 M Q whereas Davies et al. (2004) estimated a 
dynamical mass of 6.5 x 10 9 M Q within a radius of 2.5" 
(with the nucleus contributing at most ~ 15%). In addi- 
tion, G95 estimated a total stellar mass of 2.7 x 10 9 M© 
inside the ring. Thus, taking into accout the uncer- 
tainties in all the derived masses, we imposed that the 
total stellar mass fitted for the ring sections (and the 
whole ring, see below) was not greater than 3.5 x 10 9 M Q 
(~ 50% of Mdyn)- All the solutions that did not obey 
this condition were discarded. 

6 We note that a single stellar population cannot simultaneously 
fit the photometric SEDs and the spectroscopic data. Specifically, 
a 6 Myr-old population with Ay ~ 3 mag could reproduce well the 
spectroscopic data, but would not fit the photometric points. 



The derived masses, ages, and extinctions for the two 
stellar populations in the ring sections are given in Ta- 
ble [5] together with the model fits for the spectroscopic 
data (which can be compared with the observations in 
Table [5]) and bolometric luminosity. Fig. [6] shows the fits 
to the photometric SEDs. As can be seen from Tabled] 
the four ring sections have similar fitted properties. The 
young stellar population has a rather restricted age range 
(~ 5 — 6 Myr) and an average extinction of Ay ~ 10 mag. 
The fitting of the age is basically driven by the EW of Br7 
under the assumption that Ay (stars) = Ay (gas). How- 
ever, if the extinction to the gas is higher than that to the 
stars (e.g., Calzetti et al. 2000), age would be an upper 
limit. The mass in young stars for each aperture is be- 
tween ~ 250 and 500 x 10 6 M Q , comparable in most cases 
to the mass of the intermediate-age stellar population 
(see below). The only ring section with a mass in young 
stars significantly smaller than that of the intermediate- 
age population is A4 (the eastern part of the ring), as 
suggested from the lack of bright MIR emission there 
(see § 6) and the lower Br7 emission line flux (when 
compared to the other ring sections). The second stellar 
population in the ring sections has an intermediate age 
(~ 17 — 35 Myr) and is less extincted (Ay ~ 1.8 mag). 
Interestingly, these properties are similar to those of the 
majority of the 1.1 /^m-selected clusters (see § 6). 

We have also modelled the photometric SED and the 
spectroscopic properties of the whole ring, which cor- 
responds to the ring area used by G95 (see Fig. [2] 
and Table [5|). We used the Br7 flux given by G95 
(2.25 x 10~ 14 erg s' 1 cm" 2 ), and the EW of Br 7 (9.9 A) 
as measured from the mean spectrum of the three (al- 
most) non-overlapping ring sections (Al, A3 and A4) 
covered by UKIRT slits. Besides, we used the bolomet- 
ric luminosity (LrR ^Lboi — 3 x 10 11 L Q , with a 50% of 
uncertainty) of the ring derived by G95 as an additional 
datapoint to fit. The results are presented in Table [8] and 
Fig. [7] We find ages of 5 Myr and 14 Myr for the young 
and intermediate-age populations, with Ay ~ 11 mag 
and Ay ~ 2 mag, respectively. The fitted ages and ex- 
tinctions are in agreement with those fitted for the ring 
sections, not surprisingly, as the three non-overlapping 
apertures, which have all similar fitted parameters, ac- 
count for ~ 80% of the fitted stellar mass of the ring 
(Table [5]). Thus the ring sections covered by the UKIRT 
spectroscopy are representative of the ring as a whole. 
As can be seen from Table the young stellar popula- 
tion accounts for about one-third of the stellar mass of 
the ring but, interestingly, is responsible for two-thirds 
of its bolometric luminosity. 

We can independently estimate the extinction of the 
ionizing stellar population of the star-forming ring of 
NGC 7469. A number of authors (e.g., Calzetti et 
al. 2005; Alonso-Herrero et al. 2006a) have found a 
good correlation between the MIR and the extinction- 
corrected Paa luminosities for star forming regions, star- 
bursts, and LIRGs. Thus, using the 12 /im vs. Paa 
relation shown by Alonso-Herrero et al. (2006b) and 
the ring luminosity at 12.5 fim (~ 1 Jy) deduced from 
Soifer et al. (2003), we can estimate the number of ion- 
izing photons (corrected for extinction) for the ring of 
SF and compare it with the observed Br7 flux. Assum- 
ing case B recombination (Hummer & Storey 1987) and 
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TABLE 8 

Results of the x 2 -fitting of the SEDs and spectroscopic data of the ring sections and the whole ring. 
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Note. — Results of the fits using SB99 models with an instantaneous SF burst, Geneva tracks, and CAL extinction law. Br7 flux and its EW have bee: 
additional data for SED fitting. Two constrains have been imposed to the fits: extinction of the young stellar population have to be iv > 3 mag; total stellar n 
not be greater than 3.5 X 10 9 Mq (~ 50% of the dynamical mass, see text), (a) Ring sections or whole ring; (b) Parameters of the young population fit; (c) Para 
the intcrmediatc-agc population fit; (d) Age; (e) Visual extinction; (f) Stellar mass; (g), (h) and (i) Fitted Br7 flux and EW, and L^oi, respectively. These values c 
the combination of both populations; (j) Best (lowest) \ 2 value (sec Equation^; (k) Degrees of freedom, (g) and (h) fitted values can be compared with those (r 
introduced for the fit (see Table [6l. 
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Fig. 6. — Fits to the SEDs (and spectroscopic data, not shown) of the four ring sections (see locations in Fig. [2}. We imposed the 
extinction of the young stellar population to be Av > 3 mag. The observational datapoints are shown as (red) open circles with their 
corresponding 1 a uncertainties. The SED of the young stellar population is shown in blue, the intermediate-age population in orange, 
and the total (sum) fitted SED in green. The (blue) triangles, inverted (orange) triangles, and (green) squares are the model photometric 
points (SB99 SEDs convolved with ffi>T-filter throughputs) of the best fit for each population, respectively. [See the electronic edition of 
the Journal for a color version of this figure.] 



using the CAL extinction law we derive Ay ~ 13 mag 
(ABr 7 ~ 1.25 mag). This high value of the extinction con- 
firms independently (and agrees with) the existence of a 
deeply obscured young population inferred from the SED 
analysis of the ring sections and the whole ring ($ 5.1). 

5.2. Additional Spectroscopic Indicators 

The He 1 2.06 ftm/Brj line ratio has been used as an in- 
dicator of the effective temperature of the ionizing stellar 



population as well as a means of detecting the presence 
of massive stars and constraining the upper mass cutoff 
of the IMF. Under the assumptions made in this work for 
the IMF and metallicity, we can use this line ratio as an 
indicator of the age of the SF processes. For the observed 
range of the He 1 2.06 fim/Brj line ratio (0.22 — 0.29 for 
the spectroscopic apertures, see Table El and 0.29 for the 
ring, see G95) the age of the young population is approxi- 
mately ~ 5.5-6 Myr (using Fig. 3 of Rigby & Rieke 2004 
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0.2 0.5 1.0 2.0 

X (iim) 



Fig. 7. — Fit to the whole ring SED (and spectroscopic data, 
not shown). We have used the bolometric luminosity of the ring 
(G95) as an additional datapoint (see § 5.1 for details). The rest 
(constraints, symbols, and colors) are the same as in Fig|6] [See the 
electronic edition of the Journal for a color version of this figure.] 



for a Salpeter IMF, solar metallicity, instantaneous burst 
and Mup = 100 M©), in agreement with our previous 
findings (see § 5.1). However, we take these results with 
caution as Shields (1993) and Lumsden et al. (2001) have 
demonstrated that this ratio has a complex dependence 
with the effective temperature and other H II region pa- 
rameters, and cannot be used as an accurate indicator of 
the effective temperature (Rigby & Rieke 2004). 

Using SB99 models, the observed photometric CO in- 
dex (see Table is consistent with a stellar population 
with an age of 8 — 9 Myr or 50 — 100 Myr. Our results 
would be in agreement with the former interval if we 
consider that the two stellar populations fitted in § 5.1 
(~ 6 Myr and ~ 17 — 35 Myr) would generate an inter- 
mediate CO index, representative of the mixed popula- 
tions. However, as discussed by Origlia & Oliva (2000) 
the greatest uncertainty in deriving ages from CO indices 
is the stellar evolutionary tracks. This problem is allevi- 
ated if the models use observational NIR stellar libraries 
(see e.g., Forster-Schreiber et al. 2000). 

Summarizing, within the uncertainties of the models, 
the ages suggested by the other NIR spectroscopic indi- 
cators of the ring sections agree well with the ages de- 
rived from the photometric SED+Br7 fits discussed in 
the previous section, and lead to a consistent result which 
we will now place in a wider context. 

6. DISCUSSION 

6.1. Contribution of the 1.1 fim-selected Star Clusters to 
the Integrated Emission 

As we saw in § 4 and 5.1 most of the 1.1 /xm-selected 
star clusters have properties similar to the integrated 
intermediate-age stellar population inferred for the ring. 
However, the flux contribution of the clusters to the ring 
sections and the whole ring decreases with wavelength 
(Table [5]). For instance, the cluster flux contribution 
to the Al aperture is 70% (40% for the other ring sec- 
tions) in the UV and only ~ 10% in the NIR. This 
could be in part a spatial resolution effect. Since the 



spatial resolution of the UV image is lower than that 
of the 1.1 jum image, it is possible that an unresolved 
UV cluster could actually be an aggregation of clusters 
when observed at higher spatial resolution. This is the 
case for some UV and NIR clusters which have multi- 
ple cross-identifications in the higher resolution F330W 
and F606W images (see § 2.1.1). However, the UV and 
the 2.2 (im images have similar spatial resolutions, but 
the flux contribution of the clusters varies significantly 
between these two wavelengths. This suggests that dif- 
ferences in spatial resolution are not the main cause. It 
is more likely that the decreasing flux contribution (from 
UV to NIR) of the star clusters to the ring sections could 
be due to the existence of an evolved stellar population 
not spatially resolved as individual clusters, which would 
make up a significant contribution to the NIR emission, 
but almost none to the UV continuum emission. 

In terms of the stellar masses, the 1.1 /im-selected clus- 
ters only represent a minimal contribution to the total 
stellar mass (10% of the intermediate- age stellar popu- 
lation) within the ring sections. However, the images 
with the highest spatial resolutions (ACS F330W and 
WFPC2 F606W) reveal a larger number of clusters than 
those detected at 1.1 /jm (Fig.[T]and Table[3]) whose prop- 
erties are presumably similar to those described in § 4. 
Given the similar ages and extinctions of the integrated 
intermediate-age stellar population detected in the ring 
and most of the 1.1 /im star clusters, it is likely that this 
ring population is also made up of less massive clusters 
(< 1 x 10 6 Mq), not resolved in the NICMOS F110W im- 
age. As such, the majority of the 1.1 /im-selected clusters 
would be the most massive (and thus easier to detect) ex- 
amples of this intermediate-age stellar population. 

6.2. Gas/dust Distribution and its Relation with the 
Stellar Populations 

The existence of large amounts of obscuring gas and 
dust in LIRGs (and ULIRGs), and the major role played 
by the complexity of their spatial structure have been 
known for a long time (e.g., Sanders et al. 1986; Scoville 
et al. 2000; Bushouse et al. 2002; Alonso-Herrero et al. 
2006a, among many others). This is indeed the case for 
NGC 7469. The large-scale HST imaging (Malkan, Gor- 
jian & Tarn 1998; Scoville et al. 2000) of NGC 7469 shows 
that the ring has a dusty clumpy morphology embed- 
ded in a spiral-like structure with prominent dust lanes 
(Fig. U and Scoville et al. 2000; Martini et al. 2003). 
A large-scale (several kpc) stellar bar has been detected 
by Knapen et al. (2000). In addition, the CO 2-1 map 
of Davies et al. (2004) reveals the presence of molecular 
gas in the rims of the ring of SF where the dust lanes are 
and, more interestingly, forming a nuclear gas bar inside 
the ring, crossing the nucleus (see their Fig. 3). 

In moderately dusty galaxies, and specially in LIRGs, 
the spatial locations of the young stellar populations can 
be inferred from their MIR emission (see, among others, 
Soifer et al. 2001; Helou et al. 2004; Calzetti et al. 
2005; Alonso-Herrero et al. 2006b). The 12.5 /im map of 
NGC 7469 from Soifer et al. (2003, see also the Br7 map 
of G95) shows bright MIR peaks in the northeast, south 
and southwest regions of the ring (see Fig. [5]). In fact, 
the ends of the nuclear molecuah gas bar coincide with 
the position of these bright MIR peaks (obscured regions 
of ongoing SF), and these probably with the location of 
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TABLE 9 

Contribution of the 1.1 /jm-selected clusters to the ring sections and whole ring. 
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defined in § 2.2.1); (c) Fraction of the flux contributed by clusters (or ring sections) to the corresponding section (or the whole ring); (d) Fraction of the stellar 
mass (young + intermediate-age population) contributed by clusters (or ring sections) to the corresponding section (or whole ring). The number with a hash 



symbol next to a cluster label indicates the number of total additional counterparts that have been detected among all of the images (i.e., taken into account 
all wavelengths; sec Table [2] for details). 
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Fig. 8. — mpnow — "^F160W color map (equivalent to a J — H 
color map). Colors redder (darker) than mpnow — "IF160W ~ 1-03 
(typical of an old stellar population, although younger ionizing stel- 
lar populations will be bluer, see Alonso-Herrero et al. 2006a) are 
indicative of highly extincted regions. The young and intermediate- 
age 1.1 /im-selected star clusters are marked with star symbols and 
open circles, respectively. Because of the central region (r < 0'.'5) 
is strongly affected by the subtraction of the nucleus, it does not 
show any real extinction variation. [See the electronic edition of 
the Journal for a color version of this figure.] 



the Inner Lindblad Resonance(s) of the large-scale stellar 
bar (Knapen et al. 1995, Jogee, Scoville & Kenney 2005). 
The sketch shown in Fig. [10] summarizes the complexity 
of the spatial distribution of the stars, dust, and gas in 
the star-forming ring of NGC 7469. 

As can been seen from Fig. HI there is a clear spatial 
anti-correlation between the positions of the UV-optical- 
NIR clusters and the bright MIR emitting regions. The 
radio (Colina et al. 2001) and the MIR peaks in the 
ring of NGC 7469, on the other hand, are spatially co- 
incindent, as found in other LIRGs (Soifer et al. 2001). 
These displacements between the UV-optical-NIR con- 



tinuum peaks and MIR/radio bright emitting regions are 
real and correspond to distances of approximately 160 pc 
(0'.'5). 

This lack of spatial coincidence between UV-optical- 
NIR and MIR continua is in agreement with our model of 
two stellar populations to fit the integrated properties of 
the ring (and ring sections) . The UV-optical-NIR contin- 
uum (cluster and unresolved emission) traces mostly the 
mildly extincted intermediate-age (< 100 Myr) popula- 
tion, whereas the MIR and radio peaks reveal the regions 
where the most obscured and youngest (< 8 — 10 Myr) 
SF is taking place. Interestingly, this agrees with the 
fact that the 1.1 /xm-selected clusters, which mostly trace 
(and are representative of) the intermediate-age popula- 
tion of the ring (§ 6.1), are indeed located in low extinc- 
tion regions (color map, Fig. [5]) that in general do not 
coincide with the bright MIR peaks nor with the molec- 
ular gas and dust lanes in the ring (see Figs. [Hand [To]) . 
That is, most of the clusters have been detected where no 
current but rather recent (> 8 — 10 Myr) SF and low ex- 
tinction are expected, which agrees with the parameters 
found for them (~ 14 Myr and 1.5 mag of visual extinc- 
tion). The only exception are the few young 1.1 /im star 
clusters that are located in the southern part of the ring 
(where molecular gas, MIR emisssion, and NIR contin- 
uum overlap), which could be the most massive (lumi- 
nous), less-extincted examples of the integrated young 
stellar population, as already inferred from the SED fit- 
ting of the ring sections (§ 4). However, due to the small 
number of young 1.1 /im-selected clusters, they might not 
be statistically representative of the whole young popu- 
lation. 

We conclude that in the ring of NGC 7469, the UV- 
optical-NIR broad-band continuum emission reflects the 
recent SF (< 100 Myr), but not necessarily the ongoing 
(< 8 — 10 Myr) obscured SF. The location and proper- 
ties of the youngest population are only revealed by the 
MIR and radio emission, as well as the NIR Br7 hy- 
drogen recombination line. In turn, because of the spa- 
tially variable and strong obscuration in the ring of this 
galaxy, the 1.1 /zm-selected star clusters mostly trace the 
intermediate-age, low-extincted stellar population. 

7. CONCLUSIONS 
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Fig. 9. — In all the images the contours are the MIR 12.5 /im emission (Soifer et al. 2003) superimposed on the UV F218W (left, top 
panel), optical F606W (right, top panel), NIR 1.1 /an (left, bottom panel), and radio 8.4 GHz (right, bottom panel, Colina et al. 2001; 
obtained in Nov 1999) images. Note the excellent spatial correspondence between the MIR and radio emissions whereas, in general, the 
locations of brightest NIR, optical, and UV clusters (except those to the south of the nucleus, see § 6.2) do not correspond with the brightest 
peaks of MIR emission. [See the electronic edition of the Journal for a color version of this figure.] 



We have presented a multi- wavelength HST UV- 
through-NIR. imaging analysis of the circumnuclear ring 
of SF of NGC 7469, together with new ground-based NIR 
spectroscopy at four locations along the ring. These high 
spatial resolution data have allowed us to study the SF 
properties at different spatial scales, ranging from in- 
dividual star clusters (tens of parsecs) to ring sections 
(hundreds of parsecs) and the ring as a whole (diameter 
of ~ 1.6 kpc). 

At the smallest scales we have selected 30 star clusters 
at 1.1 /im. The fitting of their SEDs reveals the presence 
of two stellar populations. Most (75%) of the clusters are 
of intermediate age (8 — 20 Myr) and have low extinction 
(Ay ~ 1-25 mag). The others are younger (~ 1 — 3 Myr) 
and are moderately extincted (Ay ~ 3 mag) . Most clus- 
ters have stellar masses within the range 1 — 10 x 10 6 M , 



although more massive clusters or aggregates of clusters 
(10 - 15 x 1O 6 M ) are also found. 

On the kiloparsec scale (ring sections and whole ring) , 
the photometric SEDs and spectroscopic indicators (in 
particular Br7 luminosity and EW) were fitted with 
two distinct stellar populations: a young (~ 5 — 6 Myr), 
highly obscured (Ay — 7 — 13 mag) population, and an 
intermediate-age (~ 14 — 35 Myr), less extincted (Ay ~ 
1.8 mag) population. A single stellar population could 
not account for all the observed properties. The young 
ionizing population would be responsible for the bright 
Br7 line emission and the MIR and radio continuum 
emission peaks. This young population is generally not 
coincident with the UV-optical-NIR continuum emitting 
regions and the majority of the intermediate-age clus- 
ters. Moreover, the brightest MIR/radio regions coincide 
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Fig. 10. — Sketch of the nuclear region of NGC 7469. The approximate locations of the l.l//m-selected star clusters are marked with 
stars. The places where there is molecular gas (traced by CO J=2-l emission) are shaded in light gray, whereas the most intense MIR/radio 
emitting regions arc shaded in dark gray. The location of SN2000ft (Colina ct al. 2001, 2007) is shown as a black star. 



with the ends of the nuclear molecuar gas bar detected in 
CO. This young population of ionizing stars, although it 
only represents one-third of the stellar mass of the ring, 
accounts for a large fraction (about 2/3) of its IR lu- 
minosity. The 1 . 1 ^.m-selected clusters, although making 
up a small fraction of the ring stellar mass, are repre- 
sentative of the integrated intermediate-age stellar pop- 
ulation responsible for the UV-optical-NIR continuum 
emission. On the other hand, the younger (<~ 1 — 3 Myr) 
clusters could be the less extincted, massive examples of 
the young ionizing population traced by the MIR contin- 
uum and Br7 emission line. 

In conclusion, multi-wavelength data (UV-through- 
radio imaging and NIR spectroscopy) have allowed us to 
characterize the SF properties of both the star clusters 
and the ring of SF of NGC 7469. In the case of this LIRG, 
and other dusty galaxies, these types of observations are 
critically needed to discern and study the current (traced 
by MIR/radio continuum emission and NIR/MIR emis- 
sion lines) and the recent (detected through UV-optical- 
NIR continuum) SF taking place in the dusty environ- 
ments commonly found in LIRGs and ULIRGs. 
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